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Abstract 

We theoretically investigate the quasiparticle scattering rate F inside a vortex core in the existence of non-magnetic 
impurities distributed randomly in a superconductor. We show that the dependence of F on the magnetic field 
direction is sensitive to the sign of the pair potential. The behavior of F is quite different between an s-wave and a 
d-wave pair potential, where these are assumed to have the same amplitude anisotropy, but a sign change only for 
the d-wave one. It is suggested that measurements of the microwave surface impedance with changing applied-field 
directions would be used for the phase-sensitive identification of pairing symmetry. 
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1. Introduction 

Much attention has been focused on unconven- 
tional superconductors [1]. It is important to iden- 
tify Cooper pair symmetry for getting a clue to pair- 
ing mechanism in each superconductor. The Cooper 
pair is described by the superconducting pair poten- 
tial, which is a complex number with amplitude and 
phase. Information on its phase is crucial for identi- 
fying Cooper pair symmetry and for discriminating 
between unconventional superconductivity (e.g., d- 
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wave pair) and conventional one (s-wave pair). How- 
ever, most of experimental methods can detect only 
the amplitude of the pair potential, whereas a few 
are phase sensitive. Therefore, it is desirable to look 
for new methods capable of detecting the phase of 
the pair potential. 

In this paper, we theoretically investigate the 
quasiparticle scattering rate F inside a vortex core 
and show that it can be used as a phase sensi- 
tive probe when changing the direction of applied 
magnetic field. The field-angle dependences of the 
thermal conductivity and the specific heat have 
been discussed so far [2,3]. Instead, we will consider 
here the field-angle dependence of F inside a vor- 
tex core. The information on 7^ would be obtained 
experimentally by measurements of the microwave 
surface impedance. 
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2. Formulation 

We consider a single vortex core in a system in 
which non-magnetic impurities are distributed ran- 
domly. The quasiparticle scattering rate F is ob- 
tained from the imaginary part of the impurity self 
energy ImS. We calculate E on the basis of the qua- 
siclassical theory of superconductivity [4]. For the 
vortex bound states around a vortex core, analytic 
solutions of the quasiclassical Green's functions are 
obtained by the so-called Kramer-Pesch approxima- 
tion [5] . S is calculated from those Green's functions. 
An s-wave scattering in the Born limit is assumed 
in this paper. The pair potential is expressed as 
A(r, fcp) = A(r)(i(fcF), where the Fermi wave num- 
ber k-p denotes the position on a Fermi surface and 
the vortex center is situated at r = 0. The function 
d{kY) represents the pairing symmetry. The bulk 
amplitude of the pair potential is defined by Aq = 
|A(|r| ~> cxd)| away from a vortex core. v-p{k-p) is 
the Fermi velocity at k-p- Consider a scattering pro- 
cess of a quasiparticle from fcp to fcp on a Fermi 
surface. F is composed of an integration of the scat- 
tering process with respect to the initial state fcp 
and the final state fcp. According to procedures de- 
scribed in Ref. [6], the scattering rate F{e) for the 
vortex bound states with energy e is calculated as 



^ = |(((1 - sgn[rf(feFKfe^)]cose) 



1 \vF^{k'p)\ \d{kp) 

|sine||t;F±(fcF)| M(fc^) 



X g-«(so,feF)g-M(So,feF) 



FS' I FS 



(1) 



where 



e{kp,k'p) = e^{kp) - e,.{k'p) 



(2) 



w(s, ki 



2|rf(fcF)|Ac 

|^F_L(fcF)| 



ds' tanh 



(3) 



so(feF,feF,e) = — 



£ I \vF±{k'p 



sinel 2A2|d(A;^) 

\VF±{kF) 



2A^|d(fcF)P 



■ cos 6 



(4) 



So(fcF, k'p,e)^ so(fcF, fep, e) cos 9 

|^^F_L(feF)| 



'2Al\d{kF) 



sinO. 



(5) 



Here, F^ is the scattering rate in the normal state. 
The brackets (• • • )fs and (• • • )fs' rnean the Fermi- 
surface integrals with respect to fcp and fep, respec- 
tively, like 
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with dSp being an area element on a Fermi surface. 
We have considered a single vortex in a coordinate 
system fixed to an applied magnetic field B. clm, ^m, 
Cm are orthogonal unit vectors and cm is set parallel 
to B (cm II B). VF±{kp) is the vector component of 
vp{kp) projected onto the cm-^m plane normal to 
B. Then, ^o is defined by ^o = vp^/ttAq with [7] 
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Around the vortex, the amplitude of the pair po- 
tential was assumed to be |A(r)| ~ Aq tanh(|r|/^o) 
in the om-6m plane. O-u^kp) is the angle of VFi_{kp) 
measured from the am axis. 

The expression of F{e) [Eqs. (l)-(5)] is repre- 
sented by two coordinate systems; one is fixed to 
crystal axes and the other is fixed to the applied 
magnetic field. That is, the wave number fep is de- 
fined in a coordinate system fixed to crystal axes, 
while ^0 and the angle 6^ are defined in the coor- 
dinate system fixed to the applied magnetic field. 
Actually, there is a relation between those two co- 
ordinate systems. 

We consider a coordinate system fixed to crystal 
axes characterized by orthogonal unit vectors d, b, 
c. First, let us start with the situation where the 
two coordinate systems coincide (cm || d, Sm || b, 
Cm II c), and then rotate the vector cm (i-e., the 
direction of B) by the polar angle it/2 around the b 
axis (thus, a\i = — c and cm = d). Next, we rotate 



Cm (II B) by the azimuth angle aM around the c 
axis (see Fig. 1). In this case, we have the relation 
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In this paper, we consider a single-band super- 
conductor with an isotropic spherical Fermi surface, 
where VYik^) \\ kp and \vF{kF)\ = vf — const., 
for clarity. Using the above relation (8), one can ex- 
press VF{kp) in the coordinate system fixed to the 
magnetic field, then evaluate VF±{kF)^ and finally 
obtain the relations 
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where fcp = fcF(cos(/)fc sin^^d -f sint/ife sinfl^b -1- 
cos0kc). Thus, all quantities are represented in the 
coordinate system fixed to the crystal axes. We have 
now reached to a position where we can calculate 
Eq. (1) to investigate the dependence of 7^ on the 
field- angle au ■ 

3. Results 

In this section, we show numerical results ob- 
tained from Eq. (1). The following three types 
of Cooper pairing are considered on the isotropic 
spherical Fermi surface, (i) Line-node s-wave: 
d{kF) = I cos 20fc sin^ 6*^ I . (ii) Point- node s-wave 
[8]: d{kF) = (1 + cos4:<j)kSin^dk)/2. (iii) d-wave: 
d{kF) — cos2(/)fc sin^ 9k- Note that all these pairing 
states have gap nodes {\d{kF)\ = 0) in the 0^ — tt/A 
directions on the Fermi surface. On the other hand, 
the anti-node directions correspond to the 1/)^ — 
ones. 

In Figs. 2-4, we show the scattering rate F of 
the vortex bound states as a function of the field- 
angle au (see also Fig. 1) . Each curve corresponds to 
different energies e of the quasiparticle. In the case 
of the line-node s-wave pair (Fig. 2), F exhibits its 
minimum around the gap-node direction om — 7r/4. 
While similar behavior is seen also in the case of the 
point-node s-wave pair (Fig. 3) at the low energy 




Fig. 1. A coordinate system fixed to crystal axes character- 
ized by orthogonal unit vectors a, b, c. The applied mag- 
netic field B is rotated by the azimuth angle om in the plane 
normal to the c axis. The cm axis of the other coordinate 
system is taken parallel to B. A driving current I may be 
applied perpendicular to a plane in which B is rotated. 

e — O.IAq, a difference appears at higher energies. 
The minimum appears also around the anti-node 
direction in the case of the point-node pair, which is 
in contrast to the case of the line-node s-wave pair. 
More prominent difference appears between the line- 
node s-wave and the d-wave pair (compare Figs. 2 
and 4). 7^ always exhibits its maximum at the gap- 
node direction in the case of the d-wave pair and 
the ratio of the maximum value to the minimum one 
increases with increasing the energy e, which quite 
contrasts to the case of the line-node s-wave pair. 

These two pair potentials, (i) and (iii) above, have 
the same form in the amplitude |(i(fcF)|- The only 
difference is that the pair potential exhibits a sign 
change for the d-wave pair (iii), and not for the s- 
wave one (i). According to Eq. (1) of F, the only 
part that reflects the sign of the pair potential is 
the coherence factor [6], 1 — sgn[d{kF)d{kp)] cosG. 
Owing to this factor, arises the prominent difference 
in the field-angle dependence of F between the s- 
wave and the d-wave pair state. 

4. Discussion 

As seen in the preceding section, the field-angle 
dependence of the scattering rate F of the vortex 
bound states refiects the Cooper pairing symmetry. 
This suggests that its measurements can be consid- 
ered as a new phase-sensitive probe. F{e) of the vor- 
tex bound states is related to the flux flow resistivity 
Pi{T) ex F[e — fceT) in moderately clean systems 
[9] even when the superconducting gap has nodes 
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Fig. 2. Plot of the quasiparticle scattering rate F vs. the 
apphed-field angle om in the case of the line-node s-wave 
pair. Each curve is plotted for different quasiparticle energies 
e. The vertical axis is normalized by maximum value for each 
curve. 

on a Fermi surface [10]. The field-angie dependence 
of pi may be measured by rotating an applied mag- 
netic field within a plane normal to a driving cur- 
rent (Fig. 1) [11]. The microwave surface impedance 
measurement is a promising experimental method 
to estimate pf avoiding the vortex pinning effect at 
low temperatures [12,13] if it is possible to rotate an 
applied magnetic field there. 



5. Conclusion 

We have investigated the quasiparticle scattering 
rate F inside a vortex core under the existence of 
non-magnetic impurities. The expression of F is de- 
rived in a form where the applied magnetic field can 
be rotated in a plane. We find that the behavior 
of -T as a function of the applied-field direction re- 
flects the sign of the Cooper pairing. The difference 



Fig. 3. Plot of the quasiparticle scattering rate F vs. the 
applied-field angle ayi in the case of the point-node s-wave 
pair. Each curve is plotted for different quasiparticle energies 
e. The vertical axis is normalized by maximum value for each 
curve. 



in F between the s-wave and the d-wave pair is es- 
pecially prominent despite the same gap amplitude 
anisotropy they have. That difference arises because 
of the coherence factor which depends on the sign 
of the pair potential. Therefore, the field-angle de- 
pendence of F is sensitive to the phase of Cooper 
pair in contrast to the field-angle dependence of the 
specific heat, which is known to be phase insensi- 
tive. Elucidation of more detailed mechanism of the 
field-angle dependences in 7^ is left for future stud- 
ies. F is related to the flux flow resistivity, the field- 
angle dependence of which would be experimentally 
explored by microwave surface impedance measure- 
ments if the direction of an applied magnetic field 
can be changed there. 
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Fig. 4. Plot of the quasiparticle scattering rate F vs. the 
apphed-field angle ayi in the case of the d-wave pair. Each 
curve is plotted for different quasiparticle energies e. The 
vertical axis is normalized by minimum value for each curve. 
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